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ABSTRACT 

We explore how the behavior of galaxy cluster scaling relations are affected by flux- 
limited selection biases and intrinsic covariance among observable properties. Our 
models presume log-normal covariance between luminosity (L) and temperature (T) 
at fixed mass (M) , centered on evolving, power-law mean relations as a function of host 
halo mass. Selection can mimic evolution; the L-M and L-T relations from shallow 
X-ray flux-limited samples will deviate from mass-limited expectations at nearly all 
scales while the relations from deep surveys (10"^'* erg s~^ cm~^) become complete, 
and therefore unbiased, at masses above ~ 2 x 10^^ M©. We derive expressions for 
low-order moments of the luminosity distribution at fixed temperature, and show that 
the slope and scatter of the L-T relation observed in fiux-limited samples is sensitive 
to the assumed L-T correlation coefiicient. In addition, L-T covariance affects the 
redshift behavior of halo counts and mean luminosity in a manner that is nearly 
degenerate with intrinsic population evolution. 



Key v^rords: clusters: general — clusters: ICM 
clusters — clusters: calibration 
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^ 1 INTRODUCTION 



Encoded within observed galaxy cluster populations are im- 
portant keys to cosmic structure evolution. The decryption 
of these clues requires an accurate understanding of survey 
selection. The character of the cluster population selected 
on an observable, O, will depend strongly on how O relates 
to the underlying halo mass and how that relation varies 
over time. 

Upcoming surveys at sub-millimeter and opti- 
cal wavelengths will complem en t existing efforts in 



the X-ray d Ebeling et al.l 1 19981 : iBohrin ger et alJ 
20061: |Pierr ~ et al. 20061) and optical (jMjllcr et aj 



2002 



2005 



temperature, T — derived from observations of flux-limited 
samples. In particular, we investigate how the form and evo- 
lution of the L-T relation depend on assumptions regarding 
the covariance in these observables with mass. 

Dimensional arguments suggest power-law behaviors 
for the mean lumi nosity and te mperature as a function of 
mass and epoch (|Kaiseij Il986l ). Recent work has moved 
beyond the slopes and intercept of this relation, focus- 
ing on a stochastic model wit h log-normal scatter about 
the mean populati o n behavior llReiprich fc Bohringeil |2002| : 
IStanek et aiTl2006l : iReiprichI |2006| ). The physical origin of 
this scatter has been linked to non-adiabatic mechanisms. 



Koester et al.l |2007| ) bands. Survey cross-calibration will 



enable detailed characterizations of the inter-relationships 
among observable tracers of mass. Precisio n cosmology 
deriv e d from number counts and clustering ( Levine et al.l 
l2002l : iMaiumdar fc Mohrl |2004| : iLima fc Hul |2004| , |2005| ) 
relies on this type of interwoven tapestry. 

Although efforts to understand cluster covariance will 
ultimately be pan-chromatic, this paper focuses on the be- 
havior of two bulk X-ray properties — luminosity, L, and 
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such as cool cores (iFabian et al.l Il996 


: iMarkevitchl Il998l: 


lArnaud & Evrard 19991: O'Hara et al.l 


20061) and mergers 


(Rowley et al. 2004 iPoole et al.l 2007^. Analytic investiga- 



tions discuss the impact of entropy floors on the scatter and 
evolution of the s caling relations among luminosit y, temper- 
ature, and mass l|Boweij|l997l : iBalogh et al.ll2006l ). 

Self-similar evolution of the type envisaged by iKaiserl 

(|l986t ) has been tested with several observational sam- 
ales, including those drawn from survey s from ASCA 



pies, mcludmg those dra-w 
(|Mushotzkv fc Scharj [l997l : 



Novicki et al.l 



Newton (fVikhlinin et ahl |2002| : iKotov fc VikhlininI |2005|) 



ly s irom 

I Eooi), 



XMM- 



CHANDRA (|Lumb et aLll2004l ). and ROSAT (jFairlev et all 
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I2OO0I : iMaughan et aljlioo^ ). These studies compare the L- 
T relation among redshift sub-populations in an effort to 
find the scaling with epoch. Results range from somewhat 
stronger than self-similar scaling to slightly negative. The 
work presented here addresses a potential source of con- 
fusion for such studies. Shallow flux-limited surveys are 
mass-incomplete, while samples derived from deeper flux- 
limited samples become complete at sufficiently high masses. 
This difference in selection, essentially a redshift-dependent 
Malmquist bias, manifests effects similar to genuine popula- 
tion evolution. 

We work with two canonical flux limits, 3 x 10^^^ and 
1 X 10"^* ergs s~^ cm"'^, representative of the existing RE- 
FLEX survey and of future surveys such as ChaMP, XMM- 
LS and Spectrum X - Gamm j^ We follow the conventions 
used in IStanek et al.l (j2006l ). where L is a rest-frame, soft- 
energy band (0.1 — 2.4 keV) X-ray luminosity in units of 
lO^^ergss"^ with h = 0.7 (Ho = 100 /i^^kms^^ Mpc"^), 
and M is the mass within a sphere encompassing 200pc{z) 
in units of 10^^ Mq. Note the difference in Hubble con- 
stant conventions for luminosity and mass. We assume a 
concordance cosmology with a lower power spectrum nor- 
maliz ation that is consistent with WMAP3 (|Spergel et al.l 
[iooi): {n^^A, erg} = {0.3, 0.7, 0.8}. 

First, we discuss mass selection in the presence of scat- 
ter and flux limits and the consequences for X-ray cluster 
surveys (§2). We develop formalism in §3 for studies of cor- 
related observable properties by introducing intrinsic covari- 
ance within a joint distribution of luminosity and tempera- 
ture. Finally, we address the role of covariance in the study 
of scaling relation evolution. 



2 L-M AND MASS SELECTION 

We introduce the following notation for the logarithms of 
mass, luminosity, and temperature, respectively: fj, = InM, 
£ = InL, and t = InT. We assume that the soft-band X-ray 
luminosity of halos follows a log-normal distribution with 
dispersion, at, about a mean that varies with mass and 
epoch as 



£i5,o+PfJ. + sln[E^{z)]. 



(1) 



The normalization, £15,0, defines the present-epoch geomet- 
ric mean luminosity at the scale of 10^^ M©; p is the 
mass-scaling slope; and E{z) = ^r2m(l -I- z)^ + Qa repre- 
sents the Hubble parameter evolution for a flat metric. Our 
default model uses s = 1, appropriate for self-similar evo- 
lution in the soft X-ray band. The scatter ae is assumed to 
be mass- and redshift-independent. In what follows, we em- 
ploy values p = 1.6, ai — 0.59, and lnLi5,o= 1.81, derived 
by matching REFLEX counts to predicted cluster counts to 
the Jenkins mass function (Stanek et al. 2006). 

A sample's flux threshold, /, sets a redshift-dependent 
luminosity limit, £min{z) = \og[A-Kd\{z) f / K{z,T)], where 
the K-correction, K{z, T), is detailed in the Appendix. Tem- 
peratures are discussed in §3.1 below. For non-negligible 
scatter, ai, this sharp luminosity-selection maps to a smooth 
selection, or window, function in mass 
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Figure 1. Window functions of current (solid) and next- 
generation (solid, bold) X-ray surveys are shown along with 
the cumulative fraction, N{< z)/Ntot, (dashed) of all halos at 
IQlS/i-lM© (upper) and IQI^/i-IMq (lower). 



W{ti,z) = 



dlT{l\z,y) 



-erfc(V'min), 



(2) 



where 'P{l\z,ij) is a Gaussian distribution with mean given 
by equation ^ and constant scatter at, and ^jjmin = 
at. 

Fig. [1] displays the window functions of current and 
next-generation X-ray surveys for mass scales of the 
Coma-cluster (10^^/i~^Mq) and for halos near the lim- 
its of modern Sunyaev-Z el'dovich surveys (1O^'*^~^M0, 
iHolder fc CarlstromI [200lh . Dashed lines show the cosmic 
cumulative fraction, A'^(< z)/Ntot, of halos within redshift, 
z. For our chosen normalization ag, — 0.8, most 10^^ /i~^M0 
halos lie within z = 1, while the 10^^ /i~^M0 population 
continues growing to z = 2. The next-generation window 
function allows detection of 100% of the 10^^ h~^M.Q halos 
in the sky. The decline in the window function at z ~ 0.5 
permits detection of ~ 56% of the universal population at 
lO^'* /i~^Mq. By contrast, the shallow flux limits of current 
surveys detects only 34% and 0.06% of these totals. 

Such incomplete mass sampling leads to a Malmquist 
bias that brightens the L-M r elation relative to the intrin- 
sic (mass- limited) population (jStanek et al.|[200^ ). The log- 
mean luminosity of halos in a flux-limited sample is 



J dV{z) n{fi, z) 



£{ii,z)W{z\^j)+ate'- 



JdNf(p) 



(3) 



where dNf{^) = dV{z)n{fi, z)W{^i, z), represents the dif- 
ferential number of halos-within a comoving volume el- 
ement, dV{z)-whose luminosities satisfy the threshold, 
£min{z,T). The second term manifests the bias such that 
when ip,nin becomes negligible and W remains appreciable 
{i.e., when the median luminosity is near the flux limit), 
the geometric mean is enhanced by an amount comparable 
to the scatter, at. 

Fig- El compares the L-M relation for shallow and deep 
flux-limited surveys to the underlying, mass-limited rela- 
tion. To discern between the effects of survey selection and 
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Figure 2. Geometric mean, (L) = e^^^ , of the de-evolved luminos- 
ity for current (dotted) and next-generation (solid) flux-limited 
samples. The upper panel shows the bias relative to the underly- 
ing, mass-limited population (dashed). 



redshift evolution, we plot the de-evolved geometric mean, 
(L/i?^^), using the model value, s = 1. Samples with the 
shallow flux limit show a gradual deviation from a power 
law at the high-mass end. At extremely high masses-above 
3 X 10^^ M0-the bias disappears as the sample becomes 
mass-complete. For the deeper flux limit, surveys become 
mass-complete above 2 x 10^* Mq (see Fig. [l|, and the 
resultant L-M relation is unbiased above this mass scale. 

The deep sample probes to higher redshift than the shal- 
low sample. But simply comparing the mean luminosity at 
fixed mass between these samples, without taking into ac- 
count the evolving Malmquist bias, would produce an incor- 
rect conclusion about luminosity evolution. 



3 L-T AND CORRELATION 

Since accurate total masses are difficult to obtain, ex- 
isting studies have used X-ray tempe rature as a low- 
scatter mass proxy (|Evrard et al.l 1 19961 ). Several investi- 
gations of X-ray evolution have compared L-T relations 
am ong low- and high- r edshi ft populations. Early work 
by iMushotzky fc Schar3 l|l997h . iFairlev et all (|200(]| ). and 
iNovicki et al. (2002) find trends that are consistent with 
no evol ution ( s = ) to z < 0.8. A rece nt study by 



MaughaneTal] (|2006l ) compares the WARPS (IScharf et all 
1997|) subsample wi t h dat a from lArnaud fc Evrard ( 19991 ) 
and lVikhlinin et all (|2002l ) to probe conventional and mod- 
ified forms of self-similar evolution. They find consistency 
among the L-T slopes of differing redshift populations to 
be consistent with each other, but slightly steeper than 
the self-similar model. In a sample of 10 clusters at 2 > 
0.4, iKotov fc VikhlininI (|2005h measure evolution somewhat 
steepe r than, but consistent with, self-similarity. llTumb et al] 
l|2004l report near-self-similar evolution, and acknowledge 
the possibility of a fiux-limit bias in their eight-cluster sam- 
ple at 0.45 < z < 0.62. 

iBalogh" et al.l (|2006l ) compare observations with analytic 
models of the ICM, finding that the degree of evolution 
in the scaling relations depends strongly on that of the 




kT [keV] 

Figure 3. The X-ray luminosity— temperature relation (upper 
panels) and its dispersion it^(T) (lower panels) expected for cur- 
rent, flux-limited (3 X 10"^'^ erg cm~-^) samples and for dif- 
ferent correlation coefficients, r = (—0.7, 0.0, 0.7). Solid lines show 
the analytic expectations, points show discrete realizations from 
Hubble Volume full-sky halo realizations, and dashed lines are 
power-law fits to these discrete samples. 



L-M scatter. IStanek et al.l l|2006l ) determine the redshift- 
independent cr^|^ in the presence of a flux-limit; they also 
note that simulations suggest a weak intrinsic correlation of 
the deviations in luminosity and temperature at flxed mass. 

We next incorporate such covariance into scaling rela- 
tion analysis. After computing observables under a multi- 
variate Gaussian model, we add an explicit demonstration 
of the model usi ng halos from Hubb le Volume (HV) sky sur- 
vey catalogues l|Evrard et al.|[200^ ). Our aim is to demon- 
strate the interplay between covariance and scaling relation 
parameters for fiux-limited samples. 

3.1 Joint Luminosity- Temperature Distribution 

To calculate the geometric mean and dispersion of the L-T 
relation, we assume that the log-mean temperature of the 
intracluster medium follows virial scaling behavior. 



t ■ 



tl5 + -H' 



1 



(4) 



with normalization tis =ln(A;Ti5) and fcTis — 6.8 keV. We 
further assume that there is constant log-normal scatter 
about the mean with at = 0.10, equivalent to a 15% dis- 
persion in mass at fixed temperature. 

Using the normalized deviations. Si = — /ai and 
St = {t — i) /at, we form the joint log-normal distribution 
of luminosity and temperature at a given mass and epoch. 



where d- 



2-kR 



(5) 



(SiSj) is the correlation matrix with off- 
diagonal elements equal to the correlation coefficient, r = 
{SiSt), and R = \/l — r^. The m"^ moment of the distribu- 
tion at fixed temperature is 

^ / dVjz) J d^^ n(^, z) e-^.^/^£<-) (tj^, z) 



C^"'\t\^i,z) = 



JdNfit) 

Airv{i\z,^i,t), 



(7) 



where dNf(t) = dl/(z)d^i n(/i, z) W(m, 2)6"** The first 
moment (i.e., the geometric mean luminosity) is 



£<'>(%, 2) = W{^ji,z)l + aiRe- 



(8) 
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Figure 4. Redshift behavior of the geometric mean luminosity 
(upper), and the sky surface density (lower) of 4 keV clusters 
above a next-generation flux limit (10~^* erg s~^ cm~^). The 
solid line in all panels is our default model (s = 1, erf = 0.59, 
lnLi5_o = 1-81, r = 0). The left panels vary the correlation co- 
efficient, r = 0.7 (short-dashed), and r = —0.7 (long-dashed), 
otherwise retaining the default model parameters. The right pan- 
els hold r = and show the effects of changing the evolution and 
scatter: s = 0.95, ai = 0.72, InLis^o = 1-65 (long-dashed); and 
s = 1.1, an = 0.35, InLis o = 1-95 (short-dashed). 



where i = £ + r ae St and (pmin = (i — l) /\/2Rae. At a given 
mass, non-zero correlation effectively shifts the mean lumi- 
nosity, either enhancing or suppressing the probability that a 
halo of given mass survives the flux cut. This has important 
implications in the counts and mean masses discussed be- 
low. In the no-correlation limit (r — Q), equations Q and ((S)) 
reduce to equation ([3]), but convolved with a temperature- 
selection filter. 

From the L-T distribution's second moment, 

+e-'^-- (2RaJ+V2R^aj(Prr,,n) , (9) 

we compute the observable variance, — {l^) — (i)^. 

In Fig. O we show the expected L-T relation mean and 
dispersion for REFLEX fiux-limited clusters assuming three 
degrees of L-T correlation, r G (—0.7,0.0,0.7). Points in 
the figure represent explicit realizations of the model de- 
rived from HV sky survey samples of massive halos above 
5 X 10^^ M0. Using redshifts and masses of an all-sky 
survey, we assign luminosities and temperatures via equa- 
tions IT]) and ([4]), using two independent random Gaus- 
sian deviates, gi and g2'- gi controls the luminosity, and 
gt = r gi + R g2 sets the temperature. We then apply a 
flux limit including K-correction terms. 

Similar to the flux-limited L-M relation (Fig. [2)|, the 
L-T relation deviates from a pure power-law, although this 
deviation is extremely weak for values, r > 0.5. The scatter 
depends strongly on the correlation coefficient, with magni- 
tude varying from 0.8 to 0.4 as r varies from —0.7 to 0.7. 
Anti-correlation causes scatter orthogonal to the mean input 
relation, enhancing the magnitude of the observable disper- 
sion. Positive correlation reduces this dispersion. 



3.2 Evolution Diagnostics 

For mass-complete samples, the mean L-T relation will be- 
have as LoiT'^^/'^ pc{zY~^^'^ . For a fiux-hmited sample, the 
geometric mean soft-band luminosity at fixed temperature 
is 

fdu.n(z,a) f d£ iVUlz, ii,t) 

{eiz\t)) = ^ ^ ^ '7/ / ' ^ (10) 

where dnf{z\t) = dfin{p, z)W{ij,, z)e~^t ^'^ . The lumi- 
nosity in equation (|10|) and the sky surface density, 
{4-K)~'^dV{z) /dzdnf{z\t), are shown in Fig. U for the case 
of kT — 4 keV clusters selected with a next-generation 
(10"" erg s~^ cm"^) fiux limit. 

The left panels show models with different L-T cor- 
relation within our default model with (s, erf , InLis.o) = 
(1.0,0.59, 1.81). The survey is mass-complete at 4 keV out 
to 2 ~ 0.7, independent of r, giving counts that are identi- 
cal in this redshift range. The mean luminosities shift from 
the mass-limited mean by an amount rae{St), where (St) is 
the mean temperature deviation of the selected sample. For 
halos of fixed temperature, the steepness of the mass func- 
tion implies that (St) > 0; there are more low-mass halos 
to scatter upward into the temperature bin than high-mass 
halos to scatter downward. This effect either brightens or 
dims the mean luminosity at low redshifts, depending on 
the sign of the correlation, r. At z > 0.7, the differences in 
4 keV mass selection among the models become appreciable. 
Compared to the r = case, this selection effect drives the 
counts up and the mean masses down for the r = 0.7 model, 
and vice- versa for r = — 0.7. The magnitude of this mass shift 
is small, however, amounting to —3% and -1-6% at z = 1.5, 
respectively. 

The right panels of Fig. |4] show two zero correlation 
models tuned to produce behavior similar to the r — ±0.7 
cases. The parameter sets, (s, a^, InLis^o) = (1-10, 0.35, 1.95) 
and (0.95, 0.72, 1.65), are within the 95 % error conto u rs tha t 
resuh from the REFLEX analysis of IStanek eTeI\ (|2006l ). 
The solid line reproduces our default model, as in the left 
panels, where r = 0. The luminosity moments and sky surface 
densities of these models lie within a few percent of the 
respective covariance models shown in the left panels, but 
the driving mechanisms are different. The dimmer intercept 
and larger scatter of the s — 0.95 case force this model to 
become mass-incomplete earlier than the default case. The 
counts drop beyond z > 0.7 and the mean luminosity rises 
as the Malmquist bias is triggered. The s = 1.10 case, with a 
smaller scatter and higher intercept, remains mass-complete 
to a higher redshift, leading to larger counts and a smaller 
Malmquist bias. 

Unlike the r- = ±0.7 case, the mean halo mass selected 
by the 4 keV constraint is not shifted at high redshift in 
these models. When r = 0, halos that scatter up or down 
into a fixed temperature bin are equally likely to be culled 
by the sample flux limit. Mass estimates from weak lens- 
ing of stacked ensembles could potentially be a useful dis- 
criminatory diagnostic, but systematic effects must be con- 
trolled at the one-percent level at z~l. The upside is that 
next-generation surveys will produce large numbers of high- 
z clusters. According to Fig. HJ a 5000-square degree, deep- 
flux limit survey will yield 445 clusters with temperatures 
between 3.8 and 4.2keV, and between 0.9 and 1.1 in redshift. 
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4 CONCLUSIONS 

The multivariate properties of massive halos selected by 
X-ray flux-limited surveys are sensitive to assumptions 
about the correlation between bulk observable properties 
at flxed mass. Using a model in which the L-M relation 
is tuned to match local REFLEX counts, we show how the 
Malmquist bias in mean luminosity (which arises from mass- 
incompleteness) disappears during the transition from shal- 
low to deep flux-limited samples. Although the underlying 
model is based on power-law behavior, the geometric mean 
L-M relation from a flux-limited sample will deviate from a 
pure power-law, with a kink at the mass scale above which 
the survey becomes mass-complete. For a survey with a 
flux limit, 10~^* erg s~^ cm~^, this feature is expected to 
lie near M = 2x 10^^ h'^ Mq. 

Using a model with Gaussian covariance to describe how 
L and T jointly relate to mass and redshift, we compute 
luminosity moments and sky counts of clusters at fixed X-ray 
temperature. The slope and scatter of the L-T relation for 
bright, local samples is sensitive to the covariance between 
these properties. Placing limits on the covariance from local 
samples will require prior information on the intrinsic L and 
T variance at fixed mass, and such priors can be obtained 
from external observations or from gas dynamic simulations. 

Finally, we address attempts to extract information on 
cluster evolution from the behaviour of the L-T relation 
within a flux-limited survey. We show that non-zero L-T 
covariance affects counts and luminosity moments as a func- 
tion of redshift in a manner that is degenerate with redshift 
evolution at zero covariance. 

A deep X-ray survey of clusters, by itself, is limited to 
the information provided by counts and moments. Combin- 
ing such a sample with optical and sub-millimeter observa- 
tions offers the potential to break the evolution-covariance 
degeneracies through lensing-mass estimates and additional 
signatures that can be computed via extensions to the ap- 
proach introduced in this paper. 
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APPENDIX A: K-CORRECTION 

To convert rest frame to observed soft-band (0.5 — 2.0 keV) 
flux, we multiply the former by a correction factor obtained 
from fitting the output of an 0.3 solar metallicity mekal 
plasma model. Specifically, we use Lobs ~ K(z,kT) L-^cst 



with 

K{z,kT) = 1 + z [K-i{kT) + K2{kT) z] (Al) 

Ki{kT) = -0.209 + logfcr(1.18-0.391ogfcT) 

K2{kT) = -0.098 + log fcr(-0.092 + 0.085 log fcT), 

and kT in keV. The fit is accurate to a few percent within 
z = 2. 
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